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INTRODUCTION 

Traditionally, any attempt to describe turbulent be- 
havior in fluids starts with the Navier-Stokes equation 
(NSE) 0. However, success has at best been mixed 0- 
In this Letter, we explore the idea that NSE is not the 
unique approach to the study of turbulence, and that tur- 
bulence may be found in what we label as post-Navier- 
Stokes equations. What could modify NSE? In our opin- 
ion, the molecular nature of fluids can no longer be ig- 
nored d H H H 0. So we attempt to modify NSE 
by including quantum concepts, and in particular, in its 
simplest application, arrive at the possibility of velocity 
reversal and multi- valued velocity fields both impor- 
tant features of turbulence. Before showing these final 
results, we need to make two comments on the deriva- 
tion of NSE to justify our proposed post-NSE. 

The traditional way of deriving the Navier-Stokes 
equation (NSE) is phcnomcnological, based on the con- 
tinuum model and conservation of momentum. It may 
also be derived using the kinetic theory of structure- 
less molecules by starting with the Boltzmann transport 
equation 
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where we follow the conventional definitions from Huang 
|lfj| . For the purpose of differentiating our approach to 
arrive at post-NSE equations, we quickly comment on 
the assumptions of the derivation. 



First, collisional invariants \ are defined such that 
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These collisional invariants are x = m (mass), x — mv i 
(i = 1,2,3 momentum), \ = ^m\v — u(r, t)\ 2 (thermal 
energy), where u(r,t) — (v). 

To get the NSE, multiply the Boltzmann transport 
equation by p and integrate over all momentum, yielding 
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where p(r,t) = mn(r,t) and Pij = p((vi — Ui)(vj — Uj)). 
Our first comment is that the above equation results 
from microscopic conservation laws assuming elas- 
tic collisions of point molecules. 

To arrive at the traditional NSE, we need an explicit 
form of the pressure tensor, which is taken to be 
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(4) 



where [i is the viscosity. The above choice of the pressure 
is justified by the assumption that a fluid element, or 
microscopically, as our second comment, the particle 
of the model has no intrinsic angular momentum. 

In full component form, the Navier-Stokes equation is 
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Notice that before the introduction of the definition of 
the pressure, the conservation of momentum equation is 
exact in so far as the Boltzmann transport equation is 
valid. This equation comes only from the left-hand side 
of the Boltzmann equation. The contribution from the 
collision term disappears by virtue of conservation of mo- 
mentum and the assumption of elastic collisions. Hence, 
only the left hand side of the Boltzmann transport equa- 
tion is important to yield the conservation of momentum 
equation and NSE. We stress the two assumptions needed 
to arrive at the Navier-Stokes equation: first, elastic col- 
lisions, and second, the absence of angular momentum of 
the structureless molecules. What will happen if these 
two assumptions are no longer valid? 



CORRECTING THE NAVIER-STOKES 
EQUATION 

Suppose that each of the molecules could be found in 
any one of N states, the ground state and (N — 1) ex- 
cited states. Assume that excitations and de-excitations 
are induced by molecular collisions, which are now inelas- 
tic. Then the semi-classical analogue of the Boltzmann 
transport equation for each of N distribution functions 
will be 



where 
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where the collision term might be replaced by 
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= ^2 7mnJfm(r,p) (7) 
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7 mn is the transition probability of a particle in the 
m state jumping to the n state. J is a "jump" operator 
that will carry the conservation law the we will invoke 
in the time evolution of N distribution probabilities. K 
is a "kick" operator that allows the injection of energy 
from outside 0, 0, 0, 0] . a is the probability that 
a particle is kicked to a different momentum by external 
means. The kick operator makes the injection of energy 
into the system possible. Eq. ®is a generalization of 
our previous models [lj, [lj, [lj, . 

To calculate macroscopic averages, not only must we 
integrate over all momentum, we should also sum over 
all N states, to yield 
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We may think of © and (|10|l as the quantum correc- 
tions to the NSE. The first sum © is the driving term. 
We will call the second sum <|1(J[1 the radiative correc- 
tion because every transition is accompanied by radia- 
tion. (NSE ignores not only molecules, but also photons.) 
The radiative term represents the contribution of the in- 
ternal degrees of freedom of molecules to the macroscopic 
flow of a fluid. If all the transition probabilities are zero, 
we simply reproduce the classical NSE. This last equa- 
tion, our post-NSE, no longer assumes elastic collisions 
and the absence of angular momentum of the particles. 
The operators J and K have been defined in our earlier 
model calculations 0> they may be redefined 

with new models, but we will simplify them to show that 
even the simplest application of Equation JSJ) leads to 
novel results. 

We have had occasions to consider the radiative term 
[ill flil ITii. H3|. and will consider them later even more, 
but for now to arrive at immediate new results, we con- 
sider only the influence of quantized kicks defined by the 
operator Kf n (r,p) = f n (r,p- II) yielding 



J dp 3 pK f n (r,p) = mu + n{r 1 t)H. 
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EXAMPLES 

Remove the force F and drop the divergence of the 
pressure to obtain a non-linear equation in one dimension 
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It is the simplest application of our post-NSE but which 
remains a challenge still. 

The stationary solution is to be obtained from 
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where C is a constant and Wk denotes the k th branch of 
the Lambert W function. The function Wk(z) is a solu- 
tion of the equation we w — z in the complex plane |l5| . 
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The Lambert W function is multi-valued, making 
the stationary average velocities multi-valued. 

We will choose a toroidal geometry, and put x — 
L sin(27r6>), 6 = 0, ... 1, to ensure periodic boundary con- 
ditions. The physical model is one-dimensional, a donut 
of circumference 2-kL. One could imagine a paddlewheel 
half-stuck into the donut to provide quantum kicks to the 
fluid. If u{6 = 0) = 0, then C = Inn and we get 
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We plot Uk(0) in Figm-e^fer for k = 0, ±1, . . . ± 5. 



Real part 



Imaginary Part 




a = 0.3, m = 1, n = !, 




FIG. 1: Average velocity as a function of the angle for the 
torus geometry. The principal branch is complemented by the 
other branches which almost form a continuum. 



The solution u{x, 1) is 
u{x,t) = c 
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(19) 
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where s is defined implicitly by the second equation. 
While it is not in general possible to express u(x,t) ex- 
plicitly, we can still interpret the solution in terms of the 
characteristic curves u(x(t; s),t). 

Depending on the problem parameters a, m, II, and 
f(x), it is possible that the characteristic curves cross. If 
the characteristics cross at (x c , t), then there are multiple 
curves s\, S2, ■ ■ ■ such that x c = x(t; s\) — x(t; S2) = ■ ■ ■■ 
Moreover, the solution permits multiple values u(x c , t) = 
u(x(t] Si), t), u(x(t; S2), t), ... at the crossing point. Fig- 
ure |21 demonstrates this feature in a simple example. 





FIG. 2: The transient problem with a — 0.1, m = 1, 11 = 1, 
and f(x) = cos(x) 2 . Left: The characteristics cross at about 
t — 1. Right: The solution surface u(x,t) folds on itself and 
becomes multi- valued where the characteristics cross. 



CONCLUSIONS 



Using the method of characteristics, we now find the 
transient solutions of 



du 1 du 2 
~&b + 2~dx 



an 

m 



u(x,0) = f(x), 



(16) 



where f{x) are the initial velocity averages. Rewriting 
(|16fl as + ^ = au+ ^S, we identify its characteristic 
equations, 

x'(t)=u(t), i'(r) = l, u'(t) = au(r) + — . (17) 

m 

The boundary conditions can be parameterized as 

x (s) = s, *o(a)=0, uo(«) = /(«)• (18) 

For each fixed value of s, solving the characteristic equa- 
tions ltT7|l with initial values x(0;s) — xq(s), t(0;s) — 
to(s), u(0; s) — uq(s) yields a characteristic curve 
u(x(t; s), t(r; s)) in the solution surface u{x, i). For more 
on the method of characteristics, see for example 0. 



It appears that in our first applications of a highly 
simplified post-Navier-Stokes equation, we have arrived 
at multi- valued velocities as a function of location. They 
may well be interpreted as possible states of a turbulent 
system from which transitions to other states may be 
possible. The possibility of velocity reversal, a feature of 
turbulence, is immediately obvious. This result seems to 
be the first instance of an analytic derivation of a multi- 
valued velocity field and deserves further studies. 
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